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Abstract 

This paper reports the tenth-order QED contribution to lepton g — 2 from diagrams of three 
gauge-invariant sets VI(d), VI(g), and VI(h), which are obtained by including various fourth-order 
radiative corrections to the sixth-order g — 2 containing light-by-light-scattering subdiagrams. In 
the case of electron g—2, they consist of 492, 480, and 630 vertex Feynman diagrams, respectively. 
The results of numerical integration, including mass-dependent terms containing muon loops, are 
1.8418 (95) (a/vr) 5 for the Set VI(d), -1.5918 (65) (a/vr) 5 for the Set VI(g), and 0.1797 (40) (a/vr) 5 
for the Set VI(h), respectively. We also report the contributions to the muon g—2, which derive from 
diagrams containing an electron, muon or tau lepton loop: Their sums are —5.876 (802) (a/ir) 5 
for the Set VI(d), 5.710 (490) (a/vr) 5 for the Set VI(g), and -8.361 (232) (a/vr) 5 for the Set VI(h), 
respectively. 

PACS numbers: 13.40.Em,14.60.Cd,12.20.Ds,06.20.Jr 
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I. INTRODUCTION 



The anomalous magnetic moment g—2 of the electron is one of the most vigorously studied 
physical quantities at present, which provides a very stringent test of the validity of quantum 
electrodynamics (QED). To match the precision of the latest measurement of electron — 2 [l| 



the theory must include the QED radiative correction up to the eighth order |2h4| as well as 
the hadronic contribution (5 - 12] and the electroweak contribution 13- 15 1 within the context 
of the standard model. As a matter of fact, the largest theoretical uncertainty now comes 
from the tenth-order QED contribution which has not yet been evaluated and is given only 
a crude estimate [3] . Thus it is an urgent matter to evaluate the actual value of the tenth- 
order term. To accomplish this task we started a systematic program several years ago to 



17-221. 



evaluate the complete tenth-order contribution 

The tenth-order QED contribution to the electron g—2 consists of the mass-independent 
term and the mass-dependent terms A2 and A^ in which muon and/or tau lepton 
loop is involved, which may be expressed as 

5 



,(io) 



4 10) + Af\m e /m,) + 4 10) (m e /m r ) + A$ 0) (m e /m»,m e /m T )] (^) . (1) 



The mass-independent term A^ may be classified into six sets and further divided into 32 
gauge-invariant subsets according to the type of the closed lepton loop subdiagram. Thus far, 
;he numerical evaluation of 21 subsets has been carried out and the results were published 
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In this paper we focus our attention on the gauge-invariant set VI which consists of all 
diagrams containing a light-by-light-scattering subdiagram, one of whose photon vertex is 
external. (We call this an external light-by-light-scattering subdiagram.) Of eleven gauge- 
invariant subsets of the Set VI, eight have been evaluated previously [17|. The purpose of 
this paper is to report the evaluation of the remaining three gauge-invariant subsets: Sets 
VI(d), VI(g), and VI(h). In diagrams of Set VI(d) two virtual photon lines are attached to 
the open lepton line. This set contains 492 vertex diagrams. In diagrams of Set VI(g) one 
virtual photon line is attached to the open lepton line and the other virtual photon line is 
attached to the closed lepton loop. This set contains 480 vertex diagrams. In diagrams of 
Set VI(h) two virtual photon lines are attached to the closed lepton loop. This set contains 
630 vertex diagrams. Typical diagrams of these sets are shown in Fig. [U 
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FIG. 1: Typical diagrams of Set VI(d), Set VI(g), and Set VI(h). 

Our numerical evaluation of Feynman diagrams is based on the parametric integration 
formula 23M25| . To handle a relatively large number of diagrams systematically without 



errors, we developed an automated code-generating system called gencodeLLN that pro- 
duces FORTRAN codes for the numerical integration. It is an adaptation of the previously 
developed system for the type of diagrams without lepton loops [Tsl . Q to the diagrams 
containing an external light-by-light-scattering subdiagram. 

This paper is organized as follows. Section [III describes our scheme for numerical eval- 
uation. Section [TTT] gives the results of the numerical evaluation. Section [TV] is devoted 
to the summary and discussion. In Appendix „ we describes an algorithm for identifying 
independent set of loops on a diagram that is required for constructing the amplitude. For 
simplicity the factor (a/7r) 5 is omitted in Sees. HT1 and HTT1 



II. NUMERICAL EVALUATION SCHEME 



In this section we describe our scheme for the numerical evaluation of the diagrams of 
Sets VI(d), VI(g), and VI(h). The diagram that belongs to these sets consists of an open 
lepton line (£\) and a closed lepton loop (£2) that forms a light-by-light-scattering (l-by-l) 
subdiagram, where £\ and £2 refer to the types of leptons, i.e. electron (e), muon (m), or 
tau lepton (t). The mass- dependence of these diagrams and amplitudes is characterized by 
(£1/2) or by superscript ( ili2 \ 
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We adopt a relation derived from the Ward-Takahashi identity 

A u {p,q) ~ -q 



<?-!-0 



dp h 



dq u 

where A u (p, q) is the sum of proper vertex parts which are obtained by inserting an external 
photon vertex in the lepton lines of the self-energy function of a diagram Q in all 

possible ways. Taking account of the charge conjugation and time reversal symmetry the 
numbers of independent integrals to evaluate become 45 for Set VI(d) (see Fig. [2]), 26 for 
Set VI(g) (see Fig. ED, and 27 for Set VI(h) (see Fig. g}. 

The amplitude of a Feynman diagram is turned into an integral over the Feynman param- 
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eters assigned to the lepton and photon lines by using the parametric integral formula 
Note that the contribution of the second term on the right-hand side of Eq. fl2]) vanishes due 
to the Furry' s theorem. 

In our numerical procedure the renormalization of the amplitude is carried out by the 
subtractive renormalization. The unrenormalized amplitude Mg of a diagram Q is related 
to a finite calculable quantity AMg by appropriate subtraction terms of UV and IR di- 
vergences. These subtraction terms are prepared in the form of integrals over the same 
Feynman parameter space so that they cancel out the divergent behavior of the original 
unrenormalized integral point-by-point. 

The UV divergence arising from the 1-by-l subdiagram must be regularized, e.g., by the 
Pauli-Villars regularization. However, in Eq. ([2]) the Ward-Takahashi-summed amplitude is 
given as the differentiation of A M (p, q) with respect to q u . Therefore the divergence from the 
1-by-l loop is lifted, and the PV regularization is no longer needed. 

As a consequence the source of the UV divergences resides only in the vertex and self- 
energy sub diagrams of second and fourth order. These divergences are handled by K op- 



eration 
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26] . By definition the K operation yields the subtraction integral which is 
analytically factorizable into a product or a sum of products of lower-order quantities. It is 
symbolically denoted by 

K S Mg = L v s Y M g/s , (3) 
when S is a vertex subdiagram, and by 

K s M g = 5m^M g/s + BfMwsA , (4) 
when S is a self-energy subdiagram. Here, the superscript UV means that the leading UV 
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FIG. 2: The contribution of Set VI(d) is represented by 45 independent diagrams as listed. 
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FIG. 3: The contribution of Set VI(g) is represented by 26 independent diagrams as listed. 

divergent part is taken for the vertex renormalization constant L, the mass renormalization 
constant 5m, and the wave-function renormalization constant B, respectively. Note that 



5m^ v = 5m 2 . We also apply R subtraction [19| by which the residual part 5m = 5m — 5m 



,uv 



of the fourth-order mass renormalization constant is subtracted away to accomplish complete 
subtraction of 5m. 

Some diagrams of Set VI(d) and Set VI(g) have IR divergences. For example, the diagram 
VIg04 of Set VI (g) shown in Fig. [31 in which a photon line attached to the open lepton line 
at both ends encloses an eighth-order 1-by-l sub diagram, exhibits an IR divergence when 
;his outermost photon goes soft. These divergences are subtracted away by / subtraction 
By construction the / subtraction term factorizes as 

l S Mg = M S Lg /StU (5) 
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FIG. 4: The contribution of Set VI(h) is represented by 27 independent diagrams as listed. 

where L = L — L uv denotes the residual part of the vertex renormalization constant. 

The finite amplitude AMg obtained so far differs from the standard renormalized quan- 
tity, because the subtraction terms involve only a fraction of the renormalization constants 
relevant to the divergences. To achieve the standard on-the-mass-shell renormalization, the 
differences are collected over the diagrams of the subset, which is finite, and added to AMg. 
This step is called the residual renormalization. 

The exactly renormalized contributions of Sets VI(d), VI(g), and VI(h) to the magnetic 
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moment are given by the formulas: 

a£ 0) [VI(d) ( ^ 2) ] = AMgJ - 4ALB 2 AMjJ? + (-2ALB 4 + 5(ALB 2 ) 2 ) a g£ 



= AM%$ - 4ALB 2 - (2ALB 4 + 3(ALB 2 ) 2 ) (6) 

aj; 0) [VI(g) (£l£2) ] = AM^ 1 - 2ALB 2 AMff' - 3ALB 2 AM^ 2) + 6(ALB 2 ) 2 

= AM^ - 2ALB 2 a^ 2) - 3ALB 2 a&*> - 6(ALB 2 ) 2 (7) 
4; o) [VI(h) (£l£2) ] = AM^J - 5ALS 2 AM|; b y + (-3ALB 4 + 9(ALB 2 ) 2 ) 

= AM%$ - 5ALB 2 a^ 2) - (3ALS 4 + 6(ALS 2 ) 2 ) a^ 2) . (8) 

Here, AMvi(d), AMvi( g ), and AMyi(h) are the sum of the finite amplitude of the diagrams 
within the subsets VI(d), VI(g), and VI(h), respectively. <2 6 ll is the sixth-order anomalous 
magnetic moment containing the fourth-order 1-by-l diagram, apvb and a\y c are the eighth- 
order anomalous magnetic moments of the set of diagrams containing the external 1-by-l 
subdiagram with a virtual photon line attached to the lepton loop (IVb) or to the open 



lepton line (IVc). AMpvb arid AMpvc are their finite part defined in Ref. 261 ] . ALB 2 and 
ALS 4 are the sum of the finite part of vertex and wave-function renormalization constants 
of second and fourth order, respectively. 

The code-generating program gencodeLLN takes a one-line representation of a diagram 
as an input, and generates the numerical integration program formatted in FORTRAN. 
During this process it finds the form of the unrenormalized amplitudes, identifies the di- 
vergence structure, and constructs the UV- and/or IR-subtraction integrals. The symbolic 
manipulations concerning e.g the gamma matrix calculus and the analytic integration using 
homemade integration tables are processed with the helps of FORM [27] and Maple. 

III. RESULTS 

The numerical integration is carried out using the adaptive-iterative Monte-Carlo integra- 



tion routine VEGAS [28|. The numerical values of individual amplitudes of Set VI(d), Set 
VI (g), and Set VI (h) are listed in Tables HI HTl andlHIl respectively, for the mass-independent 
term and the mass- dependent terms in which (£i,£ 2 ) = (e,m), (m,e), and (m,t). We use 
the muon-electron mass ratio m M /m e = 206.768 282 3 (52) and the tau-muon mass ratio 
m T / m fi = 16.818 3 (27) for numerical evaluation [r) ]. 
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A. Mass-independent contribution 



Let us first consider the case in which l\ and £2 are of the same type of lepton, i.e., 
£\ — £2 — e, m,or t. This gives a mass- independent contribution to the lepton g — 2. The 
numerical values are listed in the second columns of Tabled] for Set VI(d), Table ILT1 for Set 
VI (g), and Table UTTl for Set VI(h), respectively. The values of the sixth- and eighth-order 
amplitudes and the finite renormalization constants are listed in Table IIV1 Putting these 
values into Eqs. (J6]), (J7J), and flS]), the mass-independent contributions A\ of the respective 
subsets are: 

4 10) [Set VI(d)] = 1.840 5 (95), (9) 
4 10) [Set VI(g)] = -1.591 3 (65), (10) 
4 10) [Set VI(h)] = 0.179 7 (40). (11) 



B. Mass-dependent contribution (e,m) 

The mass-dependent contribution to the electron g—2 in which the light-by-light-scattering 
subdiagram consists of the muon loop, i.e., £\ = e and £2 = m, is found from the numerical 
values listed in the third columns of Table [I] for Set VT(d), Table |TT] for Set VI(g), and 
Table UTTl for Set VI(h), respectively. The values of the mass-dependent sixth- and eighth- 
order amplitudes are listed in Table ITVl Putting these values into Eqs. dS}, ©, and ([8]), we 
obtain the mass-dependent contributions yl 2 (m e /m^) of the respective subsets: 

4 10) (m e /m M ) [Set VI(d)] = 0.001 276 (76), (12) 
4 10) (m e /m M )[Set VI(g)] = -0.000 497 (29), (13) 
4 10) (m e /m M ) [Set VI(h)] = 0.000 045 (10). (14) 



C. Mass-dependent contribution (m,e) 

Similarly, the mass-dependent contribution to the muon g—2 in which the light-by-light- 
scattering subdiagram consists of the electron loop, i.e., t\ = m and I2 = e, is found from 
the numerical values listed in the fourth column of Tables HJ [ITJ and IIHI Their contributions 
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are: 



A ( 2 W) (m^/m e ){Set VI(d)] = -7.798 (801), (15) 

4 10) (m M /m e )[Set VI(g)]= 7.346 (489), (16) 

4 10) (m M /m e )[Set VI(h)] = -8.546 (231). (17) 

D. Mass-dependent contribution (m,t) 

The mass-dependent contribution of the tau-lepton loop to the muon g— 2 is also evaluated 

and their numerical values are listed in the fifth column of Tables HJ [TTJ and IIHI The results 
are 

4 10) (m /i /m T )[Set VI(d)] = 0.081 77 (161), (18) 

4 10) (m At /m T )[Set VI(g)] = -0.044 51 (96), (19) 

4 10) (^ M / m r)[Set VI(h)] = 0.004 85 (46). (20) 



IV. SUMMARY AND DISCUSSIONS 

In this paper we evaluated the tenth-order QED corrections to the anomalous magnetic 
moments of electron and muon from the sets of diagrams, VI(d), VI(g), and VI(h). 

For the electron g—2, the total contribution is the sum of the mass-independent terms (jSJ, 
f fTUj) . and ffTTj) and the mass- dependent terms involving the muon loops f lT2|) . f fT3|) . and ffT4]) : 

/ex \ 5 

a e [Set VI(d)] = 1.841 8 (95) (-J , (21) 
a e [Set VI(g)] = -1.591 8 (65) (-J , (22) 
a e [Set VI(h)] = 0.179 7 (40) (-J . (23) 

The tau-lepton contributions to a e are more than an order of magnitude smaller than 
Eqs. (H2D, P3I>. and All]) and lie within the uncertainties of Eqs. (ED, ([22]), and ([23]). Thus 
they are negligible at present. 

For the muon g — 2, the contributions are the sums of the mass-independent terms 
ffTU]) . and f[TT]) and the mass-dependent terms involving electron loops f fT5]) . f[T6]) . and f fTT]) 
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TABLE I: Numerical values of diagrams of Set VI(d) for the mass-independent contributions (e, e) 
and the mass-dependent contributions (e, m), (m, e), and (m,t). The numerals in the parenthe- 
ses show the uncertainty of the numerical integration carried out using VEGAS. The number of 
sampling points is 10 9 times approximately 100 iterations. 
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VH01 


0.38496 


(108) 


0.0002512 


(167) 


81.810 


(140) 


0.01307 (27) 


VId02 


0.38884 


(73) 


0.0000244 


(29) 


170.126 


(91) 


0.00169 (9) 


VId03 


-0.09728 


(10) 


-0.0000264 


(4) 


4.054 


(3) 


-0.00212 (1) 


VId04 


0.31560 


(20) 


0.0000499 


(20) 


10.720 


(14) 


0.00448 (5) 


VId05 


-0.02779 


(18) 


0.0000489 


(39) 


16.903 


(20) 


0.00158 (6) 


VId06 


-0.19938 


(42) 


-0.0001046 


(46) 


-78.953 


(47) 


-0.00548 (10) 


VH07 


-0.13469 


(149) 


-0.0000488 


(66) 


-53.367 


(97) 


-0.00240 (21) 


VId08 


0.16511 


(99) 


0.0000495 


(14) 


140.628 


(93) 


0.00268 (8) 


VId09 


0.66926 


(342) 


0.0002166 


(212) 


160.196 


(211) 


0.01340 (54) 


VldlO 


-0.86487 


(129) 


-0.0003050 


(90) 


-240.503 


(138) 


-0.01865 (24) 


VIdll 


0.19659 


(187) 


0.0000178 


(131) 


82.026 


(122) 


0.00213 (33) 


VH12 


0.68550 


(148) 


0.0003757 


(439) 


104.960 


(170) 


0.02205 (61) 


VH13 


-0.84151 


(119) 


-0.0002828 


(88) 


-236.280 


(142) 


-0.01644 (22) 


VIdl4 


-0.55165 


(292) 


-0.0000066 


(105) 


-241.783 


(166) 


-0.00146 (37) 


VH15 


0.40793 


(247) 


0.0002465 


(180) 


74.465 


(151) 


0.01367 (43) 


VIdl6 


-0.86278 


(165) 


-0.0002389 


(100) 


-194.970 


(108) 


-0.01590 (31) 


VH17 


0.72320 


(103) 


0.0000980 


(53) 


335.970 


(402) 


0.00704 (14) 


VIdl8 


-0.12727 


(164) 


0.0002904 


(71) 


-105.046 


(87) 


0.01107 (22) 


VH19 


-0.68392 


(101) 


-0.0002279 


(214) 


-142.105 


(98) 


-0.01381 (36) 


VId20 


0.51681 


(295) 


0.0001113 


(72) 


293.527 


(151) 


0.00712 (33) 


VId21 


-0.84441 


(81) 


-0.0003483 


(23) 


-103.944 


(70) 


-0.02040 (8) 


VId22 


0.32927 


(138) 


0.0000666 


(71) 


37.656 


(85) 


0.00574 (20) 


VId23 


0.30074 


(175) 


0.0002004 


(126) 


90.028 


(135) 


0.01021 (31) 


VId24 


-0.06021 


(14) 


0.0000016 


(4) 


11.053 


(10) 


-0.00019 (2) 


VH25 


0.20753 


(13) 


0.0000311 


(3) 


18.408 


(11) 


0.00277 (1) 


VId26 


0.55652 


(97) 


0.0001401 


(44) 


30.649 


(32) 


0.01020 (12) 


VId27 


0.04473 


(69) 


-0.0000219 


(31) 


17.953 


(21) 


-0.00182 (9) 


VId28 


0.19223 


(194) 


0.0001718 


(146) 


78.298 


(170) 


0.00875 (33) 


VId29 


-0.55884 


(161) 


-0.0001263 


(71) 


-306.257 


(141) 


-0.00673 (21) 


VId30 


0.23422 


(35) 


0.0002785 


(65) 


29.833 


(34) 


0.01403 (11) 


VId31 


0.13852 


(41) 


0.0001240 


(86) 


50.343 


(58) 


0.00596 (12) 


VId32 


-0.65225 


(160) 


-0.0002349 


(53) 


-287.636 


(156) 


-0.01276 (20) 


VH33 


0.30989 


(145) 


0.0002517 


(96) 


38.190 


(HI) 


0.01293 (21) 


VId34 


-0.74797 


(115) 


-0.0002863 


(61) 


-158.487 


(118) 


-0.01729 (17) 


VH35 


0.36986 


(130) 


0.0002577 


(77) 


229.607 


(120) 


0.01267 (23) 


VId36 


-0.69348 


(77) 


-0.0002695 


(67) 


-81.830 


(36) 


-0.01390 (14) 


VH37 


-0.40728 


(48) 


-0.0001806 


(87) 


-80.159 


(48) 


-0.01090 (17) 


VId38 


0.59921 


(159) 


0.0001969 


(44) 


357.390 


(157) 


0.01081 (16) 


VId39 


-0.46380 


(126) 


-0.0002631 


(77) 


-78.155 


(82) 


-0.01385 (21) 


VH40 


0.11874 


(62) 


0.0000125 


(33) 


12.521 


(25) 


0.00079 (9) 


VH41 


-0.07805 


(16) 


-0.0000071 


(4) 


9.202 


(11) 


-0.00079 (2) 


VId42 


0.11040 


(83) 


-0.0001267 


(49) 


21.555 


(22) 


-0.00420 (12) 


VId43 


-0.24138 


(78) 


-0.0000597 


(39) 


18.378 


(39) 


-0.00597 (11) 


VId44 


-0.12923 


(169) 


0.0000835 


(48) 


-123.103 


(97) 


0.00486 (20) 


VH45 


0.37295 


(35) 


0.0003013 


(64) 


23.573 


(31) 


0.01616 (11) 


total 


-0.92943 (941) 


0.0007324 (757) 


37.443 (801) 


0.03080 (158) 
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TABLE II: Numerical values of diagrams of Set VI(g) for the mass-independent contributions (e, e) 
and the mass-dependent contributions (e, m), (m, e), and (m,t). The numerals in the parenthe- 
ses show the uncertainty of the numerical integration carried out using VEGAS. The number of 
sampling points is 10 9 times approximately 100 iterations. 
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-52.384 


(64) 
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(11) 


VIg02 


0.31956 


(36) 


-0.0001824 


(14) 


34.460 


(47) 


-0.00257 


(5) 


VIg03 


0.98647 


(72) 


0.0009548 


(12) 


130.839 


(34) 


0.04157 


(6) 


VIg04 


-0.29176 


(64) 


-0.0002605 


(8) 


-81.920 


(72) 


-0.01137 


(5) 


VIg05 


-2.19013 


(217) 


0.0001161 


(105) 


-185.978 


(166) 


-0.02197 


(39) 


VIg06 


0.72494 


(122) 


0.0005969 


(46) 


116.641 


(99) 


0.02417 


(17) 


VIg07 


-2.30397 


(241) 


-0.0013244 


(98) 


-178.720 


(160) 


-0.06282 


(37) 


VIg08 


1.00663 


(111) 


0.0004722 


(48) 


109.909 


(85) 


0.01977 


(16) 


VIg09 


0.77315 


(138) 


0.0000200 


(54) 


117.911 


(124) 


0.00743 


(19) 


VIglO 


-1.00933 


(221) 


-0.0004571 


(123) 


-182.581 


(183) 


-0.01757 


(43) 


Vigil 


1.84018 


(81) 


0.0013280 


(19) 


275.717 


(58) 


0.04634 


(9) 


VIgl2 


-1.22515 


(69) 


-0.0003466 


(12) 


-179.818 


(84) 


-0.02208 


(6) 


VIgl3 


2.96666 


(99) 


0.0003162 


(20) 


271.391 


(71) 


0.03606 


(10) 


VIgl4 


-1.30518 


(79) 


0.0001885 


(14) 


-172.175 


(72) 


-0.00348 


(6) 


VIgl5 


-0.60441 


(83) 


-0.0003888 


(13) 


-182.323 


(83) 


-0.01435 


(6) 


VIgl6 


1.99802 


(100) 


-0.0007980 


(22) 


275.112 


(62) 


-0.01361 


(10) 


VIgl7 


-1.85498 


(142) 


-0.0017517 


(48) 


-87.674 


(94) 


-0.07228 


(20) 


VIgl8 


0.81225 


(82) 


-0.0000606 


(21) 


58.248 


(58) 


0.00235 


(8) 


VIgl9 


0.06430 


(143) 


0.0015046 


(46) 


-88.590 


(94) 


0.05249 


(20) 


VIg20 


0.11993 


(77) 


-0.0003910 


(17) 


65.647 


(60) 


-0.01272 


(8) 


VIg21 


0.55442 


(100) 


0.0004771 


(22) 


58.433 


(83) 


0.01918 


(10) 


VIg22 


-1.87621 


(202) 


-0.0002371 


(60) 


-88.324 


(126) 


-0.02794 


(25) 


VIg23 


0.72649 


(53) 


0.0009910 


(23) 


57.310 


(59) 


0.03702 


(9) 


VIg24 


-1.42022 


(129) 


-0.0018534 


(60) 


-101.191 


(102) 


-0.06842 


(19) 


VIg25 


-0.60245 


(102) 


0.0002251 


(16) 


-164.454 


(84) 


0.00559 


(9) 


VIg26 


0.75452 


(110) 


-0.0007759 


(27) 


258.476 


(81) 


-0.02489 


(13) 


total 


1.66799 


(642) 


0.0008141 


(286) 


3.961 


(489) 


0.06914 


(93) 



and tau-lepton loop (jTBI. ( 1T9|) . and (1201: 

/a \ 5 

a M [Set VI(d)] = -5.876 (802) (-J , (24) 
<V[Set VI(g)] = 5.710 (490) (-) , (25) 
a^Set VI(h)] = -8.361 (232) (^) 5 . (26) 
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TABLE III: Numerical values of diagrams of Set VI(h) for the mass-independent contributions 
(e, e) and the mass-dependent contributions (e,m), (m,e), and (m,t). The numerals in the paren- 
theses show the uncertainty of the numerical integration carried out using VEGAS. The number 
of sampling points is 10 9 times approximately 100 iterations. 



Diagram 


AM< ee ) 


AM< em ' 




am'™' 


AAf< mt > 


VlhOl 


0.35996 (46) 


-0.0009939 


(6) 


126.906 (17) 


-0.04544 (4) 


VIh02 


0.27086 (37) 


0.0006152 


(4) 


-73.046 (10) 


0.03092 (3) 


VIh03 


-0.59951 (47) 


-0.0007719 


(4) 


50.488 (9) 


-0.03977 (3) 


VIh04 


0.26801 (45) 


0.0001968 


(4) 


-95.157 (28) 


0.01264 (3) 


VIh05 


0.13007 (35) 


0.0000478 


(6) 


1.319 (2) 


0.00311 (3) 


VIh06 


0.42117 (65) 


0.0000669 


(10) 


69.018 (48) 


0.00532 (5) 


VIh07 


3.15574 (62) 


0.0010860 


(7) 


141.662 (20) 


0.06653 (5) 


VIh08 


0.40673 (50) 


0.0000576 


(10) 


41.392 (15) 


0.00419 (5) 


VIh09 


0.02536 (48) 


0.0000732 


(8) 


17.187 (22) 


0.00327 (4) 


VlhlO 


-2.49371 (74) 


-0.0010322 


(10) 


-198.433 (40) 


-0.05948 (6) 


Vlhll 


-0.08174 (71) 


0.0001225 


(15) 


42.256 (32) 


0.00482 (8) 


VIhl2 


-1.61239 (108) 


-0.0011801 


(32) 


-94.032 (71) 


-0.06307 (16) 


VIhl3 


3.12210 (86) 


0.0000566 


(13) 


262.754 (32) 


0.01744 (10) 


VIhl4 


-3.55461 (134) 


-0.0011428 


(30) 


-209.006 (113) 


-0.07210 (14) 


VIhl5 


5.66905 (81) 


0.0021336 


(13) 


273.980 (33) 


0.12779 (10) 


VIhl6 


0.64030 (76) 


-0.0002009 


(22) 


78.458 (33) 


-0.00650 (11) 


VIhl7 


0.46696 (72) 


0.0001387 


(20) 


35.527 (46) 


0.00877 (10) 


VIhl8 


0.41308 (68) 


0.0001996 


(15) 


46.562 (34) 


0.01193 (8) 


VIhl9 


-2.53838 (76) 


-0.0008088 


(15) 


-168.962 (29) 


-0.05015 (9) 


VIh20 


1.06872 (116) 


0.0008938 


(31) 


-86.585 (72) 


0.04727 (16) 


VIh21 


2.24221 (70) 


0.0009253 


(14) 


117.262 (29) 


0.05466 (8) 


VIh22 


0.06791 (81) 


-0.0001971 


(15) 


143.881 (80) 


-0.00920 (9) 


VIh23 


-1.86591 (44) 


-0.0001763 


(4) 


-212.603 (17) 


-0.01768 (3) 


VIh24 


-0.03105 (50) 


0.0000472 


(12) 


-3.201 (4) 


0.00194 (6) 


VIh25 


-0.30141 (80) 


-0.0001359 


(21) 


51.665 (47) 


-0.00818 (12) 


VIh26 


-2.02809 (77) 


-0.0004662 


(12) 


-198.376 (59) 


-0.03184 (8) 


VIh27 


0.92582 (72) 


0.0006938 


(12) 


-98.570 (41) 


0.03794 (8) 


total 


4.54722 (393) 


0.0002482 (90) 


62.345 (231) 


0.03516 (44) 



RIKEN. T. K. thanks RIKEN for the hospitality extended to him while a part of this work 
is carried out. The numerical calculation was conducted on the RIKEN Super Combined 
Cluster (RSCC) and the RIKEN Integrated Cluster of Clusters (RICC) supercomputing 
systems. 



Appendix: Algorithm for identifying independent loops 

This Appendix describes an algorithm for identifying the fundamental set of circuits of 
a diagram following Ref. 33] which is adapted to automated handling by computers. 

In the parametric integral approach, the integrand of the amplitude is expressed in terms 
of "building blocks", B^, U, V, Aj, which are functions of Feynman parameters. B^ reflects 
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TABLE IV: Auxiliary integrals for Sets VI(d), VI(g), and VI(h). The g — 2 contribution from 
the sixth-order vertex diagram containing a light-by-light scattering subdiagram a? 1 ^ 2 ' is ana- 
lytically known for arbitrary combinations of leptons (^1,^2) 1291. l30l |. The numerical values of 
mass-dependent terms (e, m), (m, e), and (m,t) are given in [31[. The eighth-order g — 2 con- 
tributions from Group IV (b) and Group IV(c) diagrams with the lepton combinations (e, e) and 
(m, e) (ajvb> a lVb^> a iVc> ana - a iVc^) are quoted from {2, 32]. Other mass-dependent terms of these 



diagrams, ajy™\ «jvb , a iVe> ana - a iVc^' are ne wly evaluated in this paper. The forth- and second- 
order renormalization constants ALB^ and ALB 2 are related to our previous notation through 
A LB 4 = AL 4 + AB A and ALB 2 = AB 2 



Integral 



Value (error) 



Integral 



Value (error) 



(ee) 

(me) 
% 

(ee) 
°IVb 

(me) 
a IVb 

(ee) 
a IVc 

(me) 
a IVc 

ALBa 



0.371005292 . . . 
20.94792489(16) 

0.82249 (28) 
-0.41704 (375) 
-1.13891 (35) 

2.90722 (444) 

0.027930 (28) 



(em) 

(mt) 
% 

(em) 
a IVb 
(mt) 
IVb 
(em) 

(mt) 



a 



1.439445989 (77) 
0.00214283 (69) 
0.00004105 (93) 
0.006106 (31) 
-0.0001897 (63) 
-0.018233 (106) 
0.75 
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loop structure of the diagram, while Aj is related to the flow of the external momenta 



The definition of Bij is given as follows 25]. A chain diagram Q is derived from the 



diagram Q by removing all the external lines and disregarding the distinction of the type of 
lines. Suppose that a fundamental set of circuits (independent self-nonintersecting loops) of 
the chain diagram is known. Then, for % and j that label the lines of Q, is given by 

U st = }] z k ik,sik,u U = det U st , (A.l) 

*— • st 

k 

B IJ = UJ2^AU- I )st- (A.2) 

s,t 

where s and t refer to the circuits. The loop matrix ^ c takes (1, —1, 0) according to whether 
the line i is (along, against, outside of) circuit c. All lines are assumed to be appropriately 
directed. 

The circuits are found in the following way. A maximal tree T of a graph Q is a simply- 
connected set of lines that connects all vertices of Q and does not have loops. An example is 
shown in Fig. |5] where solid lines denote the maximal tree T ■ For any pair of vertices there 
is a unique path on T that links these vertices, because T is simply-connected and if there 
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FIG. 5: A choice of maximal tree T on graph Q is shown in solid lines. The chord set T* consists 
of the dashed lines. A circuit Cj is found for a line ti G T* and a path on T shown in thick lines 
that connects both ends of £j. 



were more than one path T would have a loop. T consists of n v — 1 lines where n„ is the 
number of vertices of Q. 

The chord set (or cotree) T* is the complement of T against Q, i.e., T* fl T = and 
T* U T = Q. In Fig. [5] the lines in T* are shown by dashed lines. The number of lines of 
T* is n p — n v + 1 = n\ where n p and n\ are the number of lines and loops of Q, respectively. 
Then, for each 4 G T*, there is a path V{ on the tree T which links two end-points of l± as 
shown above, li and 7^ form a closed loop Cj. These closed loops are independent with each 
other by construction, namely, li ^ Cj if i ^ j. They form a fundamental set of circuits. 

For the diagrams without closed lepton loops (q-type diagrams) as described in Ref. 18], 
the fundamental set of circuits are found rather trivially. In this case the maximal tree T is 
chosen as the set of lines that forms the open lepton line, and thus the chord set T* consists 
of the photon lines. Therefore, each circuit is identified by a photon line and a string of 
lepton lines that connects both ends of the photon line. 

For a general diagram, a maximal tree of the diagram can be found in the following way. 
Choose a vertex t> of Q, and starting from t> , extend the tree by adding a line adjacent 
to the vertex of the tree, of which the other vertex is not yet included, until all vertices of 
Q are visited. During this process a unique path V(vo,v a ), i.e. a string of lines in T that 
runs from Vq to v a , is also found for each vertex v a . A path on T connecting an arbitrary 
pair of vertices v a and Vb is then found by joining the paths V(vb,v ) and V(v ,v a ) where 
'PivbjVo) = 7 ? ~ 1 (fo, is the string of lines of V(vo,Vb) in reversed order. Note that the 
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duplicated lines in the paths are to be eliminated. 
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